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Antiferromagnetic interlayer coupling in Fe/c-SiFe/Fe sandwiches and multilayers 
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We report the first ab-initio study of the interlayer exchange coupling in Fe/c-FeSi/Fe sandwiches and Fe/c- 
FeSi multilayers. We perform several structural studies, which show the stability of the CsCl arrangement 
seen experimentally for the spacer. We find antiferromagnetic coupling between the Fe slabs for spacer 
thicknesses smaller than about 15 A, for both sandwiches and multilayers structures. We also study the 
effect of pinholes, interface roughness or structural misconfigurations of the spacer on the sign and magnitude 
of the exchange constant J. We finally show that the asymptotic behavior of J is determined by a flat band 
of the c-FeSi spacer, located at the M point in the Brillouin zone. 
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Magnetic multilayer devices (MMD) can be roughly 
classified according to the conducting nature of the 
spacer material, which leads to markedly different be- 
havior of the exchange coupling constant between the 
magnetic layers, J, as a function of the thickness of the 
spacer, z [El. The first class comprises those systems 
whose spacer is a metal. In this case, J(z) shows os- 
cillatory behavior, whose period is typically of a few 
Angstroms; superimposed to these oscillations, J(z) also 
decays as 1/z 2 , becoming negligible after several tens of 
Angstrom. The second class is composed of those de- 
vices where the spacer is a semiconductor. Now, J(z) is 
frequently antiferromagnetic (AF) and its magnitude de- 
creases exponentially with a decay length of at most two 
or three Angstrom. (Fe/c-FeSi) MMD stand out among 
all such structures as having a metallic spacer in the 
brink of becoming a semiconductor. Consequently, (Fe/c- 
FeSi) MMD might display the characteristic features of 
crossover, or even critical, behavior between two differ- 
ent physical regimes. Indeed, the exchange coupling con- 
stant of these devices is always AF, but has a fairly large 
decay length, becoming negligible at spacer thicknesses 
more proper of a metallic than of a semiconducting spacer 
|§-j7|. The details of growth and the crystal structure of 
(Fe/c-FeSi) MMD, in addition to the unique behavior of 
J(z), might also make these systems strongly attractive 
for manufacturers interested in the design of spin polar- 
ized transport devices. Indeed, the iron silicide spacer, 
which possesses the CsCl structure (c-FeSi) |,|,|, can be 
formed by inter-diffusion of Fc and Si slabs of appropriate 
thicknesses grown epitaxially [0J|]. LEED and AES ex- 
periments, which have been applied successfully to study 
epitaxial Fe/Si/Fe (001), have shown that the perpendic- 
ular interlayer distance in this bcc-like structure remains 
constant at wl.43 A, very close to the values for pure 
bulk bcc Fe §. 

Despite the large amount of experimental information, 
and the plausible technological relevance of these MMD, 



ties have been performed up to now. Only the asymp- 
totic region of exponentially decaying antiferromagnetic 
exchange has been qualitatively analyzed in an article 
by Vries et al. |(| using both Bruno's ideas JD and a 
model by Shi and coworkers [^0| . We report in this Let- 
ter a thorough study of the structural stability, the ex- 
change constant J, and the magnetic moment distribu- 
tion of Fe/c-FeSi/Fe(001) trilayer and Fe/c-FeSi multi- 
layer devices as a function of the spacer thickness z. We 
show how interstitial defects such as pinholes or inter- 
face roughness affect the sign and magnitude of J for thin 
spacers, by studying several mixed interfaces. We have 
extended Bruno's theory Q to account for the effect of 
different effective masses on the exchange constant. This 
has allowed us to make an excellent fit to the asymptotic 
behavior of J, which shows that the relevant contribution 
to the exchange constant comes from a band in the c-FeSi 
spacer located at the M point of the Brillouin zone, and 
not at the X point, as speculated previously |rjJio|. 

We have determined the spin-polarized electronic 
structure using a scalar-relativistic version of the fc-space 
TB-LMTO method |ll|] developed in the atomic spheres 
approximation. We have used two different versions of 
the generalized gradient approximation (GGA) for the 
exchange and correlation potential: the Perdew-Wang 
(PW) H and the Langreth-Mehl-Hu (LMH) Q, and 
benchmark them against preliminary results obtained in 
the Local Spin Density Approximation (LSDA) jl4|] . We 
model the system according to the existing LEED struc- 
tural analysis of MBE-grown Fe/SiFe/Fe (001) sand- 
wiches briefly mentioned above 0,^] . We have performed 
our calculations using both the experimental interlayer 
distance (1.43 A) and the lattice parameters which come 
from a total energy minimization of Fe and c-FeSi bulk 
materials, and obtained overall qualitative agreement be- 
tween both sets of calculations. For the sake of clarity, we 
discuss here the results obtained using the experimental 
distance, unless otherwise stated. We employ the super- 



Fe7/c-Fe(„_i)Si n /Fe7 in the (001) direction, with n = 1 - 
6. Two successive trilayer sandwiches are then separated 
by enough layers of atomic empty spheres to have each 
individual trilayer decoupled from the rest. We have also 
investigated the effect of the thickness of the Fe substrate 
on the stability of the interlayer exchange coupling, and 
found that the results are not altered for Fe slabs thicker 
than the selected case of seven layers. The electronic 
and magnetic structures are calculated using an increas- 
ing number of fc-points in the irreducible Brillouin zone. 
Convergence is obtained for 90 fc-points, results for 132 
fc-points being quantitatively the same. 
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FIG. 1. Total energy difference J = Ef - Eaf as func- 
tion of the number n of Si atomic monolayers in the spacer 
for Fe7/c— Fe(„_i)Si n /Fe7(001) multilayer devices. The in- 
set shows J for trilayers. Open squares denote values of J 
obtained using the theoretical lattice constant and the LMH 
functional. 

We have first made structural analyses to test that the 
experimental result of inter-diffusion of a thin Si film into 
the Fe slabs, that creates the c-FeSi structure, can be un- 
derstood in terms of the energetics of the different plausi- 
ble atomic arrangements in the spacer. We have therefore 
computed the total energy for Fe7/Si3Fe2/Fe7 TD with 
different layer configurations in the spacer, which mod- 
ify the c-CsCl structure. We have checked, for instance, 
that the Fe7/Si/Fe/Si 2 /Fe/Fe7 arrangement is 300 mRy 
higher in energy than the Fe7/Si/Fe/Si/Fe/Si/Fe7 con- 
figuration, and that in the former case the F alignment 
is slightly more stable than the AF one. 

We have further performed molecular dynamics sim- 
ulations, relaxing the atomic positions in the z direc- 
tion, for Fe 9 /c-(FeiSi 2 )/Fe9 and Fe7/c-(Fe 2 Si 3 )/Fe7 TD 
using the SIESTA code. SIESTA is a Density Func- 
tional Theory code based in pseudopotentials which al- 
lows structure optimizations of the whole system which 
are beyond the capabilities of the all-electron TB-LMTO 
method, particularly when surfaces are present. Besides, 



tures composed of iron atoms in different environments 
fill . We have used a minimal basis set and the PW 
functional for exchange-correlation jl^. Pseudopoten- 
tials were generated with the Troullier-Martins method 
Jl?| , with 4s 1 3d 7 and 3s 2 3p 2 valence configurations for 
Fe and Si, respectively. We have found that the forces 
in the theoretical bcc-like configuration just make the 
interstitial Fe atoms move slightly into the FeSi spacer, 
and expand a bit the Si-Si distances, the structure being 
very stable otherwise. The obtained exchange constants 
J are also similar within 10 % to those obtained with the 
TB-LMTO method in the unrelaxed structure, therefore 
lending further support to our theoretical study. 

We pass on now to discuss the behavior of J(z), defined 
as the difference between the total energies of F and AF 
alignments of the Fe slabs, for thin TD with perfect inter- 
faces. Fig. 1 shows the exchange constant as a function of 
the number n of Si atomic monolayers, which is a measure 
of the thickness of the spacer. We always find a positive 
J of the order of some mRy, which decreases with z in a 
non-monotonic fashion, has a bump for n = 3 and van- 
ishes for n w 5 - 6. The bump might correlate with a small 
protuberance which is seen experimentally at exactly the 
same spacer thickness ||. It is not an artifact produced 
by the TB-LMTO code because we also find it when we 
use SIESTA. Its height depends, on the other hand, on 
the choice of the interlayer distances, becoming smaller 
when we use theoretical or relaxed values for them (see 
the inset in Fig. 1). The experimental peak seen for n=4 
seems to be somewhat sample dependent [ fl8[ , and we find 
no traces of it in our simulations. The exchange constant 
becomes negligible for spacer thicknesses larger than five 
or six atomic Si layers, which correspond to a distance of 
about 15 A. Because the error bars coming from our cal- 
culations are of the order of the values obtained, we can 
only honestly conclude that J does not have any large os- 
cillation for such thicknesses. Experiments actually show 
that for z larger than «13A, J(z) follows the asymptotic 
behavior of a semiconductor, with a large decay length 
of about 3.6A This overall theoretical behavior for thin 
trilayer devices, that was also obtained within the LSDA 
is therefore in nice qualitative agreement with the 
experimental results by Vries et al. 0] . 

We find that the exchange constant obtained for mul- 
tilayers follows the same trend as that found for tri- 
layer structures, its magnitude being roughly a factor 
of two larger (see inset of Fig 1). This fact can be 
qualitatively explained in terms of simple Heisenberg-like 
physics, making use of an analogy with spin chains: the 
energy to flip a spin in a molecule composed of two atoms 
is half that required to do so in an infinite chain. We also 
find that for given z, J slightly oscillates as a function of 
the thickness of the iron slabs, as seen in Co/Cu MMD 
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FIG. 2. Distribution of magnetic moments for 
Fe7/c-Fe3Si4/Fe7 (001) sandwiches in the (a) ferromagnetic 
and (b) antiferromagnetic configurations. Black bars corre- 
spond to the magnetic moments of Fe whereas white bars 
correspond to Si moments. 

Vries and coworkers observe F-type contributions su- 
perimposed on top of the AF-like behavior in their Kerr 
hysteresis loops for TD with spacer thicknesses smaller 
than about 6-7 A || . We present now results for several 
TD where we simulate a variety of atomic misconfigura- 
tions inside the spacer or at its interfaces with the iron 
slabs, because we wish to test whether they can induce 
such ferromagnetic coupling. The first one consists of the 
Fe7/Si/Fe/Si2/Fe/Fe7 TD which was introduced a few 
paragraphs before. It is supposed in this case that the 
Si layers have not diffused completely into the Fe slabs 
to form the c-FeSi structure. We find that J is ferromag- 
netic, with J=-3.49 mRy. We present next the case of 
dense arrays of thin pinholes or rough interfaces for TD 
with n = 1, 2, and 3 (e.g.: up to about 8 A). We simulate 
them by doubling the cross section of our unit cell, which 
now contains two atoms per atomic plane. For n=l, we 
exchange one atom in the Si layer by another in the last 
layer of one of the iron slabs. We create this way a pe- 
riodic array of thin pinholes which leads to a strongly 
ferromagnetic J of -73 mRy (about thirty times larger in 
magnitude than the corresponding one for a pure c-FeSi 
spacer). We subsequently look at two different configu- 
rations for n=2. In the first one, the Fe and Si atoms are 
arranged again in such a way that there are thin bridges 
along the spacer as before. The exchange constant is 
then also negative, with a value of -33 mRy. We make 
the second configuration such that one of the Si layers 
has no iron atoms, which means that one of the inter- 
faces is rough while the other one is perfect. We find in 
this largely asymmetric case that J is equal to -0.7 mRy, 
which is of the same order of magnitude than its antifer- 
romagnetic counterpart. We finally present for n = 3 the 
case where both interfaces are rough while the central Si 



is again largely ferromagnetic, J = -5.1 mRy. We there- 
fore conclude that atomic misconfigurations gcncrically 
give rise to very strong ferromagnetic couplings. A com- 
parison of our results with those presented by Vries |(| 
lets us infere some further conclusions about the struc- 
ture of these sandwiches, namely: (a) it is very likely that 
for TD with thicknesses smaller than about 3 A there is 
a process of diffusion of iron into the Si layer (or vice 
versa), so that the Si layer is disrupted; (b) for TD with 
spacer thicknesses larger than about 4 A, diffusion takes 
place but leads to the formation of a c-FeSi spacer with 
interfaces of high quality. 

We come to comment now on the profiles that we find 
for the magnetic moments, which we show in Fig. 2 for 
TD with n=4. First, we have that the c-FeSi spacer dis- 
plays tiny magnetic moments of order 0.05 /is at the 
interface with the iron slabs, with whom they couple an- 
tiferromagnetically. Second, the absolute values of these 
magnetic moments are almost identical for both align- 
ments. Third, iron atoms in the vicinity of the inter- 
face with the spacer have reduced magnetic moments 
(~ 1.46/is) due to the hybridization of their orbitals 
with those belonging to c-FeSi atoms. Surface effects, on 
the other hand, induce enhanced magnetic moments at 
the Fe external layers, (M~ 2.97/^b). Finally, magnetic 
moments at the center of the Fe slabs slightly oscillate 
around the bulk value (~ 2.20/is). We find similar be- 
havior for multilayer devices, apart from the obvious fact 
that no surface effects exist for them. 
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FIG. 3. (a) Details of the band structure of bulk c-FeSi. 
The band drawn with a dashed line has a band edge at the M 
point and a large fiat section around the X point; (b) Total 
density of states of bulk c-FeSi (in states/eV), showing the 
band edge at the Fermi level and a sharp peak located about 
0.3 eV above it. 

The last part of this letter is devoted to study the 
asymptotic behavior of J(z). We have extended Bruno's 
theory to account for the different effective masses of elec- 
trons located at the relevant bands of the spacer and of 
the majority and minority populations of the iron slabs, 
a.nd used it to studv the contribution to ,T(z) comine from 



authors []6|,g_0| have speculated that the main contribu- 
tion to J(z) should come from a band centered at the X 
point, which gives rise to a sharp peak in the density of 
states, see Fig. 3. We have generated the band struc- 
ture of bulk ferromagnetic bcc iron and c-FeSi using the 
TB-LMTO code and used them to fit such a band, as 
well as those in the iron slabs closest to the Fermi en- 
ergy, to parabolae. This has allowed us to extract band 
edges Eg and effective masses m* , which we have plugged 
into our extension of Bruno's formula for J(z). We have 
found that the exchange constant decays in an oscilla- 
tory fashion. This behavior, which is due to the negative 
effective mass of the relevant band at the spacer, can not 
be easily reconciled with the experimental data. We have 
therefore turned our attention to the contribution of the 
same band at the M point, which now gives rise to a 
three-dimensional band edge right at the Fermi energy. 
The naive estimate for decay length A would in this case 
be infinity. If such an hypothesis was confirmed, the ex- 
change constant would decay as a power law, displaying 
therefore the typical behavior of critical phenomena jl9| . 
We nevertheless find that additional contributions to A 
come from the fact that the E^'s of the bands at the iron 
slabs are not placed at the Fermi energy, a fact which 
gives rise to finite effective barriers (we only find power 
law decay when one of these edges is finely tuned to zero 
energy). We have fitted again the £ 9 's and m*'s of the 
relevant bands at the M point and introduced them into 
the formula for J(z). We have extracted this way a decay 
length of 3.6 A and more generically an overall behav- 
ior in excellent agreement with the experimental results 
available to us. Fig. 4 shows how J (z) fits the exponential 
e -z/3.6^ rpj^g i a yout of this figure purposely mimics that 
of Fig. 5 in Ref. || with which should be compared. We 
wish to stress that all parameters entering into the ex- 
pression for J(z) have been obtained independently from 
a band structure calculation, which has left us no room 
for fitting. 
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FIG. 4. Asymptotic behavior of the exchange constant 
as a function of the thickness of the spacer. Solid line is the 
result from our amended version of Bruno's theory; dashed 
line is a fit to p~ 2 / 3 - 6 _ The inset shows the logarithm of J 



In conclusion, we have performed a thorough study 
of the structural stability and magnetic properties of 
Fe/c-FeSi magnetic trilayer and multilayer devices, us- 
ing a combination of density functional methods. We 
have studied how dense arrays of pinholes and interface 
roughness, or other structural misconfigurations affect 
the AFM character of the exchange constant. We have 
deduced that the process of diffusion which is used to 
grow these MMD is very effective and leads to interfaces 
of very high quality. We have made use of a slight mod- 
ification of the theory by Bruno to show how it is the M 
point which governs the behavior of J(z) for large thick- 
nesses of the c-FeSi spacer. 
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